Generalized tomographic maps 
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We introduce several possible generalizations of tomography for quadratic surfaces. We analyze 
different types of elliptic, hyperbolic and hybrid tomograms. In all cases it is possible to consistently 
define the inverse tomographic map. We find two different ways of introducing tomographic sections. 
The first method operates by deformations of the standard Radon transform. The second method 
proceeds by shifting a given quadric pattern. The most general tomographic transformation can 
be defined in terms of marginals over surfaces generated by deformations of complete families of 
hyperplanes or quadrics. We discuss practical and conceptual perspectives and possible applications. 
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I. INTRODUCTION 

Most of classical applications of tomography are based 
on light propagation along optic rays (implicitly assumed 
to be straight lines). Standard Radon transform theory 
guarantees that a measurement of the absorption of light 
beams travelling in dielectric media in straight Hues al- 
lows the complete reconstruction of the matter density 
of these media. Indeed, the original Radon transform [l| 
maps functions of two variables in the plane onto func- 
tions of one real variable on a line and one variable on 
a circle. The crucial property is that the transform is 
invertible and continuous [3, Q . 

There exist several generalizations of the Radon trans- 
form. See, e.g., i and 0. Further generalizations can 
be motivated by physical observations: for instance, if 
the function on the plane is a probability density, its 
Radon component is a family of probability densities of 
one random variable on the line, parameterized by a vari- 
able living on a circle A tomographic approach in a 
similar framework was applied to a free classical particle 
moving on a circle Q, where the phase space is a two 
dimensional cylinder. 

In quantum mechanics the Radon transform of the 
Wigncr function Q was considered in the tomographic 
approach to the study of quantum states [lo| and ex- 
perimentally realized with different particles and in di- 
verse situations fill . [T3| . Other experiments have 
been proposed [14| and the whole field is in continuing 
evolution, also in view of its relevance in genuine quan- 
tum mechanical problems and quantum information re- 
lated topics. Good reviews on recent tomographic ap- 
plications can be found in Ref. [la, with emphasis on 
maximum likelihood estimations [l6|, that enable one to 
extract the maximum reliable information from the avail- 



able data. 

A further development, extending the analysis to in- 
corporate more general symplectic transforms, was pre- 
sented in (TtI and the mathematical mechanism at the 
basis of the mapping of true density states onto tomo- 
graphic probabilities was elucidated in [l8| . There is an 
interesting relation between the Radon map of Wigner 
functions and the formalism of star product quantiza- 
tion [l^ : symplectic tomograms are indeed the Radon 
components of the Wigner function and this enables one 
to define a procedure aimed at determining the marginal 
probability densities along straight lines in phase space. 
The knowledge of all these marginals makes possible the 
reconstruction of the Wigner function in the quantum 
case and of the probability density in the classical case. 

The generalization of tomographic maps to curved sur- 
faces opens new perspectives in the applications of to- 
mography both to quantum and classical systems. Some 
attempts to study marginals along curves other than 
straight lines were introduced in Ref. . Very recently, 
optical "accelerating" Airy beams were observed j2l| : 
these beams could be used to perform a tomographic 
map over parabolas in phase space. A generalization of 
tomography to this kind of applications requires a gen- 
eralization of the Radon transform. 

The aim of this article is to study generalizations of the 
Radon transform to multidimensional phase spaces and 
to marginals along curves or surfaces. Most of the gener- 
alizations of the Radon transform proceed by considering 
geodesic submanifolds of a given Ricmannian manifold. 
We develop here a different approach, that can be applied 
to the Radon components of the probability densities of 
classical particles in phase space, and construct the cor- 
responding tomographic maps. 

This article is organized as follows. In Sec. |TT]we re- 
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FIG. 1: Tomography on the plane; {q,p) G R^,d G R, 6* £ § 
(unit circle). 

view the standard tomographic application of the Radon 
transform on the plane. In Sec. Illll we consider the gen- 
eralization to arbitrary dimensions. A deformation of the 
Radon transform with applications to elliptic and hyper- 
bohc problems is presented in Sec. IIVI In Sec. |V] we 
introduce a new type of transform involving hyperbolic, 
elliptic and parabolic quadrics. The transform is defined 
by translations of a basic pattern. Finally in Sec. IVII 
we discuss the relevance of our results for future applica- 
tions. 



II. TOMOGRAPHY ON THE PLANE 

Let us consider a function f{q,p) on the phase space 
{q,p) G of a particle moving on the line g G K. The 
Radon transform, in its original formulation, solves the 
following problem: reconstruct a function of two vari- 
ables, say f{p, g), from its integrals over arbitrary lines. 

In the {q,p) plane, a line is given by the equation 

X^^iq-lyp = 0, (1) 

with (/i, t^) ^ (0,0). Thus, the family of lines has the 
manifold structure K. x §, with § the unit circle, d = 
AT/VAi^n^ G M and pi/iy = tan6', 6* G S (sec Fig. [J). 

It is possible to write the Radon transform in affine 
language (tomographic map) [l|, 0] 

ujf{X,fi,u) = {S{X-^q-vp)) 

= I f{q,p)S{X - ^iq~iyp)dqdp, (2) 



where 6 is the Dirac function and the parameters 
X,^i,v € M. 

The inverse transform of ^ reads 

fil.p) = / c.,(^,M,-)e'(-----)^^^. (3) 

The homogeneity of ujf{X^ /i, v) 

ujf {XX, A/i, Xi^) = w/ (A, (4) 

VA G M, A 7^ 0, is a direct consequence of ([2]). If the 
function f{q,p) is a probability density function on the 
phase space of a classical particle, i.e. 

/(9,p)>0, / fiq,p)dqdp=l, (5) 

the function w/(A, fi, v) is also nonnegative and is called 
symplectic tomogram or Radon transform of the distri- 
bution function f{q,p) (in analogy to the Fourier trans- 
form of a function). The Radon transform contains the 
same information on the state of the particle evolving on 
the phase space as the initial distribution function. In 
summary, the tomograms 

ujf{X,^,iy) >0, / ujf{X,fi,iy)dX = 1, (6) 

form a family of density functions that depends on the 
two real parameters ii and v. 

HI. TOMOGRAMS ON HYPERPLANES 

The above construction can be generalized to higher 
dimensional spaces in a straightforward way. Let us con- 
sider a function /(g) on the n-dimensional space q G M". 
Is it possible to reconstruct the function / from its inte- 
grals over arbitrary (n — l)-dimensional linear subman- 
ifolds? The answer to this question is positive and pro- 
vides a generalization of the original Radon transform. 

A generic hyperplane is given by the equation 

X-n-q^Q, (7) 

with A G M and /i G R". Due to homogeneity, this 
family of hyperplanes is an n-dimensional manifold dif- 
feomorphic to R x S"~^, because any hyperplane can be 
characterized by its unit normal vector and its dis- 

tance to the origin |X|/|^|. Note that this manifold is not 
diffeomorphic to M" because the sphere S"~^ is compact. 
The Radon transform is given by 

LVf{X,^l) = (5(A-M.g)) 

= / /(g)<5(A - M • g)d>. (8) 

When 71 = 2 Eq. ^ is recovered. 
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FIG. 2: Diffeomorphism of the plane: G ^ (x,y) 



The inverse transform of 



reads 



(27r)" 



The homogeneity of u; / {X, ji) 



(9) 



(10) 



VA S M, A 7^ 0, is a direct consequence of ([8|). If the 
function f{q) is a probabiUty density function on R" 



/('Z)>0, / /(g)d"g = l, (11) 
also the tomograms ujf{X,fi) are probability densities 
cjf{X,n)>0, [ ujf{X,fi)dX ^1, V^ieM" (12) 



and the family of tomograms depends on the n real pa- 
rameters fi. In quantum mechanics such construction was 
applied to Wigner functions providing a center of mass 
tomography |22l |. 



IV. TOMOGRAMS ON HYPERSURFACES 

A simple mechanism that allows nonlinear generaliza- 
tions of the Radon transform is the combination of the 
standard transform with a diffeomorphism of the under- 
lying R" space. Let us consider a function f{q) on the 
n-dimensional space q S R". The problem is to recon- 
struct / from its integrals over an ri-parameter family of 
submanifolds of codimension one. 

We can construct such a family by diffeomorphic de- 
formations of the hyperplanes (in the a; G R" space) 



X — fi ■ X = 0, 



(13) 



with X G R and jj. G R". Let us consider a diffeomor- 
phism of R" 



q e 



X = ip{q) G R". 



(14) 



The hyperplanes (jl3p are deformed by ip into a family of 
submanifolds (in the q space) 



X - ^l ■ ifiq) = 0. 



(15) 



The case n = 2 is displayed in Fig. [2j {q,p) G R^ ^ 
{x,y) ^ ^{q,p) GR"^. 

Given a probability density f{x) on the x space, the 
Radon transform can be rewritten as 



ujf{X,fi) 



where 



{S{X - fi • x)) 

[ f{x)6{X - fi ■ x)d''x 

j{v{qmX - M . ^(g)) J(g)d"g, (16) 



(17) 



dxi 






dqj 




dqj 



is the Jacobian of the transformation. 
Observe now that 

f{x)d^x = f{^{q))J{q)d"q, 



(18) 



whence f{q) = f{ip{q))J{q) is a probability density. 
Therefore the tomograms are given by 

u;f{X,ti) = {6{X - ■ ^{q))) 

f{q)S{X - ■ ^{q))d"q, (19) 



with X G R and /X G R". 

The inverse transform follows by (O: 

f{q) = hv{q))J{q) 



c.H^,/.)J(<z)e'^(^--''»^, (20) 



(27r)" 



dqj 



i(X-li-ip{q)) 



with a modified kernel 
X(g;X,/x) = J(g)e'(^-''-'^(«» 

Therefore, a probability density distribution on R" 

/(<Z)>0, / /(g)d> = l, (22) 



(21) 



produces tomograms u!f(X,n) that are probability den- 
sities 



ujfiX,n)>0, ujf{X,ii)dX = 1, V^iGR". (23) 

The family of tomograms depends on the n real param- 
eters /X. We can now consider different applications of 
these deformed generalizations of the Radon transform. 
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A. Circles in the plane 



B. Hyperbolas in the plane 



In the punctured {x, y) plane without the origin (0, 0), 
the conformal inversion 
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P 



{x,y) = v{q,p) = 1 2 , 2' 2 , 2 

maps the family of lines 

X — — vy = Q 
into a family of circles 



centered at 



C 



(JL M 

\2X' 2XJ 



(24) 

(25) 
(26) 
(27) 



and passing through the origin (see Fig. [3]). When X ~ 
they degenerate into lines through the origin. 
The Jacobian reads 



J{q,p) = 



d{x,y) 



d{q,p) 



1 



(28) 



whence the transformation is a diffeomorphism of the 
punctured plane of functions / G L^(M^). The singular- 
ity of the transformation at the origin (0,0) is irrelevant 
for tomographic integral transformations because it only 
affects a zero measure set. 
Equations (Hn])-® become 



ojfiX,^,,.)^{^5l^X 
= [ f{q,p)d(x 



vp 



q- 



vp 



q2 _^ p2 



dqdp (29) 



In the (x, y) plane the family of lines 
X — px ~ vy — 
is mapped into a family of hyperbolas 

X-^-vp^O, 

q 

with asymptotes 
by the transformation 



X 

P = — , 

V 



{x,y) = ip{q,p) = ( -,p 



(31) 
(32) 

(33) 
(34) 



For fi > the hyperbolas are in the second and fourth 
quadrants, while for /x < they are in the first and third 
quadrants (see Fig. |4]). When /x = or ^/ = they 
degenerate into horizontal or vertical lines, respectively. 
The Jacobian reads 



Jiq,p) 



d{x,y) 



d{q,p) 



1 



(35) 



whence the transformation is a diffeomorphism in the cut 
plane without the axis (0, y). 
Equations ([Tg)) - (pn)) become 



ujf{X,ii,v) = [S[X---vp 



[ f{q,p)s(x-^-up)dqdp{36) 

\ q J 
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and 



The tomograms (|T9|) , ([29]) and ([36]) have the homogeneity 
property ([TO]) . 



C. Hyperboloids in R" 

The generahzation to higher dimensions of tomo- 
graphic maps that can be given in terms of quadratic 
expressions is straightforward. Let us consider for exam- 
ple the following tomographic map 



6{X-n-q-,,{q,p))f{q,p)d^qdy, 



where p and q are vectors in R" and 



Vjqjpj. 



(38) 



(39) 



This map corresponds to a deformation of the standard 
multidimensional Radon transform by means of the fol- 
lowing diffeomorphism of M^"\ ljj {((7,p) : qj = 0} 

whose Jacobian is 

d{x,y) 



(40) 



d{q,p) 



(41) 



The inverse map is given by 

dXdJ'pLdJ'v 



fiq.p) = 



2n + l 



ujf{X,fi,iy) 



(2^)2^ 



This corresponds to the higher dimensional generaliza- 
tion of the Bertrand-Bertrand tomography Q . 

Note that when n = 2, by interchanging the role of X 
and —/X, one recovers the same distribution of hyperbolas 
in the plane analyzed in the previous subsection IIV Bl 

Although the above generalizations might be very use- 
ful for light rays tomograms, all of them involve inte- 
gration over unbounded submanifolds. One would like to 
generalize the Radon transform to marginals defined over 
compact submanifolds. that are bounded on a compact 
domain around (p, q) . This case will be investigated in 
the following section. 



V. TOMOGRAMS ON QUADRICS 

Let us now look for a different generalization of to- 
mograms. We shall consider marginals along compact 
quadrics. This can be achieved by shifting and scaling a 
given quadric pattern 



X = {q- p,B{q - p)), 



(43) 



where J5 is a non-degenerate symmetric operator with 
respect to the scalar product {x,y) = x • y. A new gen- 
eralization of the tomographic map can be defined by 

uO}{X,y^-B)^ I f{q)S{X-{q-p,B{q-p)))d^q. 



(44) 

Eq. (|44p defines a completely different type of transform, 
supported on the quadrics defined by equation (|43p . It is 
easy to show that the inverse map is defined by 



m 



xe' 



j(X-(g-M,S(<?-A'))) 



(45) 



Indeed, by applying the definition of tomographic map 

mi 



IdetBl 



dX d''fi S^-i1-^'M'l-^^)))^^(^x, ^; B) 



IdetSl 



dXd^fie 



i(X-(g-M,-B(g-M))) 



X J d"e<5(X-(e-/i,B(^-/i)))/(0, 

which after integration over X yields 
IdetSl 



(46) 



/(O 



IdetBl 



'[(e,S«)-(<7,S9)+2(9-e.BM)] 



/(Oc'[«^^«)-(«-^')l<5" (q - = fiq)- (47) 



The meaning of the above tomographic map depends on 
the features of B. If we assume that B is strictly posi- 
tive (elliptic case) , this map corresponds to averages of / 
along the ellipsoids defined by Eq. ([43]) . In particular if 
all the eigenvalues of B are equal to it corresponds to 
integration over spheres centered at fi of (squared) radius 
X/lP', namely. 



h^{q~pf^X. {X>Q) 



(48) 



Note that, in the two dimensional case, the distribu- 
tion of circles is different from that obtained by the trans- 
form defined by diffeomorphisms in Sec. IIVI There, the 
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FIG. 5: Tomography on circles of center /x and (squared) 
radius X/b^ on the plane. 



family of tomograms were defined only on circles pass- 
ing through the origin, including straight lines (circles 
of infinite radius). Here, we are taking into account all 
possible circles of finite radius in the plane (see Fig. [5]). 
This corresponds to trajectories of particles moving in 
a plane under the action of a constant magnetic field. 
From a practical perspective, this new tomographic map 
would make possible a different practical implementation 
of tomography. 

When B has both positive and negative eigenvalues 
this correspond to hyperbolic tomography with averages 
of / along the hyperboloids defined by Eq. (|43l) . e.g.. 



(49) 



In the case of degenerated B forms we have to consider 
a hybrid transform. B can then be decomposed into a 
non-degenerated bilinear form and a linear form. In this 
case the tomography of the linear components should be 
treated as the standard Radon transform, whereas the 
non-degenerate variables should transform as above. Let 
us consider for example a simple three-dimensional case 
with 



B 




(50) 



In this case we can define the following tomographic map 
UJs{X,il-B)= [ d\f{q) 

x5{X-{qi- fiif - (92 - ^i2f " [i-iqz) , (51) 
with inverse transform 

dXd^[i 



fil) = 



27r3 



-iOfiX,^i;B) 



xe 



:(X-((7i-/^i)^-(g2-A'2)^-/J3g3) ^52) 



VI. CONCLUSIONS AND PERSPECTIVES 

Let us discuss the main findings of this article, both 
from mathematical and physical perspectives. From 
a mathematical viewpoint, the generalizations of the 
Radon transform introduced here enable one to consider 
marginals defined over submanifolds that are not nec- 
essarily geodesic submanifolds in Riemannian spaces or 
Lagrangian submanifolds of symplectic manifolds. These 
transforms define tomograms over compact submanifolds 
and can be more suitable for physical applications, be- 
cause the practical implementation of a tomogram can 
only achieved in local terms. In this new framework the 
recovery of a local value (of a probability distribution on 
phase space in the classical theory, or of a Wigner dis- 
tribution in the quantum case) only involves integration 
over manifolds that do not reach infinity. In a previous 
article 0], we considered the tomography of a classical 
particle moving on a circle, which required the defini- 
tion of marginals over the helices on a cylinder. Now, in 
the light of the new transforms just introduced, we have 
the possibility of performing tomography over compact 
submanifolds even for classical systems that evolve in un- 
bounded domains. This is a significant conceptual step 
forward. 

Physically, the new reconstruction formulas enable one 
to generalize the measurement procedures of the mat- 
ter density of an object. In a material medium with a 
strongly inhomogeneous refractive index, the radiation 
beams (light beams, sonic beams or matter waves) would 
propagate along curved lines and yet yield complete in- 
formation on the matter distribution by means of gen- 
eralized Radon transforms. For illustrative purposes our 
examples focus on two-dimensional situations [see e.g. 
Eqs. (|34p and (|^^ ] but the approach wc propose is 
more general and valid in R". 

In quantum optics the new "non-linear" Radon trans- 
forms can be easily extended to the quantum domain by 
using the Weyl- Wigner map. This will be discussed in a 
future publication. The results of this article show that 
the reconstruction of the Wigner function using optical 
or symplectic tomography based on straight-line Radon 
transform can be extended to situations in which the 
marginals in phase space are measured for curved hy- 
perbolas or ellipses. In particular, parabolic tomography 
could be implemented with the recently observed accel- 
erated Airy beams [2l|. 

Novel physical applications of tomography have at- 
tracted increasing attention during the last few years. 
Recent applications involving neutrinos, e.g. to get a 
mapping of the Earth inner density [23| . do not require 
new concepts of tomographic maps. However, neutrino 
tomography of gamma ray bursts and massive stellar col- 
lapses 2J] might require generalized tomography. In par- 
ticular, gamma rays tomography that made possible the 
discovery of asphericity in supernovae explosions |25j or 
imaging of astrophysical sources [26j can involve non- 
linear trajectories of gamma rays due to strong gravi- 
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tational Icnsing effects. In those cases, generalizations of 
tomographic maps hke the ones considered in this paper 
are necessary. 

Acknowledgments 

We thank J.F. Carinena and F. Ventrigha for helpful 
discussions. V.I.M. was partially supported by Italian 
INFN and thanks the Physics Department of the Uni- 



versity of Naples for the kind hospitality. P.F. and S.P. 
acknowledge the financial support of the European Union 
through the Integrated Project EuroSQIP. The work of 
M.A. and G.M. was partially supported by a cooperation 
grant INFN-CICYT. M.A. was also partially supported 
by the Spanish CICYT grant FPA2006-2315 and DGIID- 
DGA (grant2006-E24/2). The work of V.I.M. was par- 
tially supported by the Russian Foundation for Basic Re- 
search under Project No. 07-02-00598. 



[1] J. Radon, Uber die bestimmung von funktionen durch 
ihre integralwerte Idngs dewisse mannigfaltigkeiten, Bre- 
ichte Sachsische Akademie der Wissenschaften, Leipzig, 
Mathematische-Physikalische Klasse, 69 S. 262 (1917). 

[2] F. John, Plane waves and spherical means: Applied to 
Partial Differential Equations (Wiley Interscience, New 
York, 1955). 

[3] R. S. Strichartz, American Mathematical Monthly 89, 
377 (1982). 

[4] S. Helgason, Ann. of Math. 98, 451 (1973); Groups and 
Geometric Analysis (Academic Press, Orlando, 1984); 
The Radon Transform (Birkhauser, Boston, 1980). 

[5] I. M. Gel'fand and G. E. Shilov, Generalized Func- 
tions: Properties and Operations, Vol. 5 (Academic 
Press, 1966). 

[6] O. V. Man'ko and V. I. Man'ko, J. Russ. Laser Res. 18, 
407 (1997). 

[7] M. Asorey, P. Facchi, V. I. Man'ko, G. Marmo, S. Pas- 
cazio and E. C. G. Sudarshan, Phys. Rev. A 76, 012117 
(2007). 

[8] E. P. Wigner, Phys. Rev. 40, 749 (1932); J Moyal, 
Proc. Camb. Phil. Soc. 45, 99 (1949); M. Hillary, R. F. 
O'Connell, M. O. Scully, and E. Wigner, Phys. Rep. 106, 
121 (1984). 

[9] J. Bertrand and P. Bertrand, Found. Phys. 17, 397 
(1987). 

[10] K. Vogel and H. Risken, Phys. Rev. A 40, 2847 (1989). 
[11] D. T. Smithey, M. Beck, M. G. Raymer, and A. Faridani, 

Phys. Rev. Lett. 70, 1244 (1993). 
[12] G. Zambra, A. Andreoni, M. Bondani, M. Gramegna, M. 

Genovese, G. Brida, A. Rossi, and M. G. A. Paris, Phys. 

Rev. Lett. 95, 063602 (2005); M. Genovese, G. Brida, M. 

Gramegna, M. Bondani, G. Zambra, A. Andreoni, A.R. 

Rossi, M.G.A. Paris, Laser Physics 16, 385 (2006); G. 



Brida, M. Genovese, F. Piacentini, Matteo G. A. Paris, 
Optics Letters, 31, 3508 (2006). 
[13] C. Kurtsiefer, T. Pfau, and J. Mlynek, Nature 386, 150 
(1997). 

[14] G. Badurek, et al, Phys. Rev. A 73, 032110 (2006). 
[15] Quantum State Estimation, Lecture Notes in Physics Vol. 

649, edited by M. G. A. Paris and J. Rehacek (Springer, 

Berlin, 2004). 

[16] P. Banas, J. Rehacek, and Z. Hradil, Phys. Rev. A 
74, 014101 (2006); Z. Hradil, D. Mogilevtsev, and J. 
Rehacek, Phys. Rev. Lett. 96, 230401 (2006). 

[17] S. Mancini, V. 1. Man'ko and P. Tombesi, Quantum Semi- 
class. Opt. 7, 615 ( 1995). 

[18] V. I. Man'ko, G. Marmo, A. Simoni and F. Ventriglia, 
Open Sys. & Information Dyn. 13, 239 (2006). 

[19] O. V. Man'ko, V. I. Man'ko and G. Marmo, Phys. Scr. 
62, 446 (2000). 

[20] M. A. Man'ko, V. I. Man'ko and R. V. Mendes, J. Phys. 

A 34, 8321 (2001). 
[21] G. A. Siviloglou, J. Broky, A. Dogariu and D. N. 

Christodoulides, Phys. Rev. Lett. 99, 213901 (2007). 
[22] A. S. Arkhipov and V. 1. Man'ko, Phys. Rev. A 71, 

012101 (2005). 

[23] C. Kuo, H.J. Crawford, R. Jeanloz, B. Romanowicz, G. 
Shapiro and M. L. Stevenson, Earth and Planetary Sci- 
ence Letters 133 (1995) 95. 

[24] S. Razzaque, P. Meszaros and E. Waxman, Phys. Rev. D 
68, 3001 (2003). 

[25] K. Maeda et al. Science Express (online edition of Sci- 
ence), 31 January 2008. 

[26] S. N. Zhang, G. J. Fishman, B. A. Harmon and W. S. 
Paciesas, Nature 366 245 (1993). 



